Switching the sign of Josephson current through Aharonov-Bohm interferometry 
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We investigate the DC Josephson effect in a superconductor-normal metal-superconductor junc- 
tion where the normal region consists of a ballistic ring. We show that a fully controllable tt -junction 
can be realized through the electro-magnetostatic Aharonov-Bohm effect in the ring. The sign and 
the magnitude of the supercurrent can be tuned by varying the magnetic flux and the gate voltage 
applied to one arm, around suitable values. The implementation in a realistic set-up is discussed. 
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The broad interest in mesoscopic physics has recently 
spurred a renewed impulse in the context of Josephson 
effect. Because of its large spectrum of applications to 
nanotechnology, the art of manipulating the supercur- 
rent is presently under the spotlight 1] . Josephson field- 
effect transistors, for instance, have been proposed Q 
and realized both with semiconductors and carbon nan- 
otubes Q. A growing interest is nowadays devoted to 
the issue of supercurrent sign reversal: a tt-J unctions 
state, i.e., a Josephson current flowing in the direction 
opposite to the phase difference between the supercon- 
ductors, has already been obtained with ferromagnet- 
superconductor (FS) junctions 4]. In these systems the 
sign of the current flow depends on the F-layer thickness, 
which cannot be varied during an experiment. In view 
of large-scale technological applications, the tunability of 
a system plays instead a crucial role, and the realization 
of controllable 7r-junctions represents a major challenge. 
To this end, two approaches have been explored so far: 
tunable 7r-states have been realized either by driving far 
from equilibrium the junction quasi-particle distribution 
function through current injection from additional termi- 
nals [5] , or by exploiting the electron-electron correlation 
in a quantum dot, where the electron number can be 
tuned through a gate Q . 

In this letter, we present a different operational princi- 
ple to implement a controllable 7r-junction, which is based 
on the electro-magnetostatic Aharonov-Bohm (AB) ef- 
fect in a ballistic normal (N) ring coupled to two super- 
conducting (S) leads. The proposed set-up is shown in 
Fig. [TJ the ring is threated by a magnetic flux , and 
a gate voltage Vg is applied to one of the arms. Exploit- 
ing the magnetic or the electrostatic AB effect, one can 
tune the transmission of the ring and therefore manip- 
ulate the magnitude of the Josephson current. Here we 
show that, if both the magnetic and the electric field are 
suitably applied, also the sign of the supercurrent can be 
changed. This leads to a fully controllable system op- 
erating at equilibrium, which exhibits potential for the 
implementation of low-dissipation transistors, quantum 
interference devices, as well as reduced- noise qubits 
Ideal candidates for the realization of the proposed set- 




FIG. 1: Scheme of the Aharonov-Bohm Josephson junction. 
Light gray region indicates the normal metal ring, whereas 
dark gray denotes the superconducting leads. 

up are intermediate/long SNS junctions fabricated with 
InajGai-a; (with x > 0.75) ;8:] or InAs [§| rings, since 
they allow ballistic electron transport and lack of a Schot- 
tky barrier when contacted to metals, e.g. superconduct- 
ing Nb. 

We describe the system by setting the superconducting 
order parameter as A{x) ~ |A|e^'-^/^ in the left lead, 
A{x) = |A|e+'^/2 in the right lead, and A(a;) = in 
the N region, where x is the phase difference between 
the S electrodes. The Josephson current J{x) of an SNS 
junction with a N-region characterized by a Fermi veloc- 
ity Vp and a length L can be written [l^j El as a sum 
of two contributions Jd(x) and Jc(x)) arising from the 
discrete Andreev levels and the continuous spectrum re- 
spectively. Exphcitly, Jd,c(x) = (2eEL/h)jd,c{x), where 
El = hv^ /L is the energy associated with L (here the 
length of each ring arm), and 

MX) - Et-M^)(-|f) (1) 

Jcix) = —T- / d£ In 2 cosh — . 2 

In Eqs. energies are expressed in terms of E^, 

so that |Al| = \A\/El, e = {E - Ef)/El is the en- 
ergy variation with respect to the Fermi level Ep of leads 
and N-region, and /3l = El/UbT is the inverse temper- 
ature. The Andreev levels ep(x) S [0; |Ai|] in ^ are 
determined by the solutions of D{e\x) ~ 0, where the 
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function 

Die; x) = det[l - A{s)r\{x)S{e)rA{xW{-e)] (3) 

accounts both for the dynamics in the N-region, through 
the particle (hole) scattering matrix S(e) (S*(— e)), and 
for the Andreev reflections at the S contacts through 
the matrix rA{x) — diag(e*-^/^, e~*-^/^) and the am- 
plitude A{e), which equals exp(— 2i arccos(£/|Ai|)) for 
|e| < I Ail and {e - ^ - / {e + ^ for 

jej > |Al|. Finally, the function H appearing in ^ reads 



i7(e;x)=axlni?(e;x). 



(4) 



Although general, Eqs. ([I])-(l2l) do not allow a straight- 
forward evaluation of J(x) for an arbitrary S-matrix, 
for the determination of the Andreev levels may be a 
formidable task. In the limit of short junction |Al| <C 1, 
however, Beenakker showed that the contribution ([2|) 
vanishes, and that the evaluation of the Andreev levels 
in ll]) is considerably simplified. He thus proved [l^ that 
in this limit and for a symmetric S-matrix, the Joseph- 
son current acquires a simple expression in terms of the 
transmission coefficients T^'s of S. In a typical experi- 
mental realization of the set-up of Fig. [l] however, the 
short junction regime |Ai| ^ 1 is not achieved; as ob- 
served above, one rather has |Ai| > 1. Furthermore, 
in the presence of a magnetic fiux, the ring S-matrix is 
not symmetric in general. It is therefore important to 
analyze the behavior of the Josephson current also be- 
yond the short-j unction limit, and without assuming the 
symmetry of the S-matrix. To this purpose, it is suit- 
able to rewrite Eqs. (HI)-© in a different way. Here we 
briefly sketch the main steps of our strategy, based on 
analytic continuation of Eq. ^ in the complex plane. 
We first observe that, due to causality [13], the scatter- 
ing matrix S(e) has an analytic continuation S(z) in the 
upper complex half-plane Im(z) > 0. Then, elementary 
properties of holomorphic functions lead to prove that 
Hk{z; x) = IuDrXz; x) and Hx{z] x) = \nDA{z; x) 
with 

D^{z; x) = det [1 - A{z) r\{x)S{z)r a{x){S{- z*))*] (5) 
D^{z; x) = det [1 - A{z) r^(x)(S(z*)) Vl(x)S(-z)] (6) 

are analytic continuations of fl]) for Im(z) > and 
Im(z) < respectively. Here A{z) is the continuation 
of A{e) with a branch-cut in z e [-|Al|; |Al|]. The 
relation iJR(e + iO+) = H^{£ -f iO") for |e| < |Al| 
stemming from ([5]) -([6)), together with equality d^Ep = 
—d^D{ep;x)/deD{ep;x) for the Andreev levels, allows 
to rewrite Eq. ^ as 
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(X) = ^J2J2[ tanh(^)i?,(z,x)rfz, (7) 

T> A ,.1 i 11 
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where (^^) is a small semicircular contour around Sp, 
in the upper (lower) half-plane. Similarly, the relation 



ifR(e-|-iO+) = i7X(e + «0~) for \e\ > \Al\ stemming from 
©-(HI) leads to cast Eq. ([2]) in the same form as Eq. (O, 
where the sum over F^ (r^) is replaced by a contour 
along the upper (lower) real axis range e g [|Ai|;oo]. 
Finally, the analyticity of iJR(z;x) and Ha{z;x) allows 
to merge the above two contributions into 
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^ Reln(DR(zw;„;x)) , (8) 

m=0 



where w„i — (2m + l)7r//9i are the Matsub- 
ara frequencies in units of El/H. The relation 
Hfi^{~iwn;x) ~ H^{iwn]x) f^as been used. 

Equation ([8]), combined with Eq. ([5]), gives the Joseph- 
son current in terms of the N-region scattering matrix S. 
Although equivalent to the original Eqs. HI)-®, Eq. ^ 
does not require to determine the Andreev levels, yielding 
a major simplification in computing J(x) for non-short 
junctions [13.]. We shall now show the rich physical sce- 
nario arising from a non symmetric S-matrix. In order to 
illustrate this, we focus here on the case of a single con- 
duction channel, where the S-matrix is a 2 x 2 unitary 
matrix. Any U(2) matrix can be univocaly written as 



S(£) = 



p{e) -iris) 
-iT*{e) p*{e) 



(9) 



where ]p(£)]^ + ]r(e)]^ = 1 and ^ is a radius (i.e., ]^] = 1). 
Inserting Eq. ^ into Eq. ([3]) one obtains 

D{e; x) = 1 + A'(e)e2'*(^) - 2A(e)e'*(^) x (10) 
[Re {p{e)p*{-e)) + \T{e)T{-e)\ cos(x - Xo(e))] , 



Notably, Eq. (|10p shows that the phase difference x is 
renormalized by a shift xo, related to the phase of the 
amphtude T{e) — ]r(e)] exp{i(j)r{£)) through the relation 



Xo(e) = (j)rie) + 4'ri~e). 



(11) 



From Eqs. (jlip it follows that xo is non vanishing 
only if r is a complex number. At first sight, this 
property might seem a quite general feature of any 
quantum interferometer. This is, however, not the case, 
for the entries Sy of the scattering matrix obey the 



micro- reversibility relation 1^ ^ijls = Sji]_B, where 
B is the magnetic field. Thus, for any interferometer 
operating in the absence of magnetic field, the S-matrix 
is symmetric, i.e., r(e) is real (see Eq. ([9])). One thus has 
Xo = and the system always behaves as a 0-junction. 
By contrast, ii B ^ 0, the time-reversal symmetry (TRS) 
is broken, opening the possibility to the asymmetry 
of the S-matrix, and hence to a xo 7^ 0- Importantly, 
-B 7^ is a necessary, but not sufficient condition to 
have a 7r-state: in a symmetric ring threaded by a 
uniform B, for instance, S is still symmetric and one 
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has a 0-junction. If, however, the interferometer is 
suitably designed, an appropriate tuning of its parame- 
ters can lead to Xo — inducing a switching to a 7r-state. 

The electro-magnetostatic AB interferometer (see 
FiglD) is an illuminating example to describe this ef- 
fect. The explicit expression for the AB ring S-matrix 
can be obtained with standard techniques |14| by com- 
bining the scattering matrices describing the y-junctions 
with the propagation matrix along the ballistic arms. For 
simplicity we neglect here band curvature, fringe field 
effects, and spin-orbit interaction. Each contact is de- 
scribed by the 3x3 F-junction matrix [iBl, with 



entries Sj^ = — cos 7, 



■-'13 



-'21 



^31 



5^3 = (cos7- l)/2, and S^g = S^2 = (cos7 + l)/2, 



->22 

where the parameter 7 S [0; 7r/2] accounts for the contact 
transmission, with 7 = 7r/2 describing a fully transmit- 
ting contact and 7^0 the tunnel limit ■ The trans- 
mission of the two contacts will be assumed equal. The 
propagation along the two arms leads to the AB interfer- 
ometry effect; right movers, for instance, acquire a phase 
exp{i{kL — <j)/2)) along arm 1, and exp{i{kL + u + <p/2)) 
along arm 2, where (j) = ^Tre^s/h and u = eVa/EL (see 
FiglD). Similarly for left-movers. After lengthy but stan- 
dard algebra, one can compute the S-matrix ((9]) and the 
phase shift xoi which turns out to fulfill 



ReTp(e)Rerp(— e) — Imrp(e)Imrp(- 
\Tp{e)Tp{-e)\ 



cos xo(e) 

\Tpyt)Tp\-t.)\ 

(12) 

where iTp{e) = X]s=± exp[is{e + k'p)], C± — cos[(u± 
(/))/2], and Kp = kpL + u/2, with hp the Fermi wave vec- 
tor. Notice that xo is independent of the contact trans- 
mission. For the purely electrostatic AB effect (</> = and 
u ^ 0), we obtain xo = 0, as expected from the preser- 
vation of TRS. Moreover, Eq. (fT2|) yields xo = also 
for the purely magnetic AB effect {u — and 0^0): 
this is due to the additional relation Si2|s = S12I-B, 
which holds in a symmetric ring if u = 0. Although this 
relation breaks down if the ring is realized asymmetri- 
cally, the arm length is not a tunable quantity. However, 
Eq. ([T^ indicates that a much simpler way to achieve a 
controllable 7r-state is to combine the magnetic flux and 
the gate voltage: when both 7^ and u 7^ a phase 
Xo 7^ arises. Notice that, in this case, a Josephson 
current can flow even if x = 0; if TRS is broken, the 
amplitude of processes bringing Cooper pairs from right 
to left lead arc not necessarily compensated by those re- 
lated to the opposite direction. 



By substituting the S-matrix into Eq. (jlOp , and by per- 
forming the analytic continuation Z?r as discussed above, 
the Josephson current is determined by Eq. ([8]). We shall 
now discuss the results. Figure [2] displays the switching 
of J(x) from to 7r-junction behavior, for a set-up with 
|A| = 2El and contact transmission 1/2 (7 = 7r/3) which 
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FIG. 2: (color online) Sign switching of the Josephson current 
J for the set-up of Fig. [T]with |A| = 2El, kpL = 60 and 
contact transmission 1/2 (7 = 7r/3). a) j(x) for = and 
u = (solid), for <f) = Q and u — O.Stt (dotted), and for 
</> — 0.8n and u — O.Stt (dashed); b) j{x ~ 0) vs. <j!> for 
different values of the gate voltage: u = 7r/20 (dotted), u = 
tt/2 (dashed), and u = 1.2ti (solid); c) j{x = 0) vs. u for 
different values of magnetic flux: cj> = Q (solid), cj> = O.Itt 
(dotted) and (j> — O.Qtt (dashed). 



corresponds, e.g., to a InGaAs ring with L ~ 600 nm [8| 
contacted to Nb leads, at zero temperature 17|. In par- 
ticular. Fig. mja) shows J as a function of the phase dif- 
ference X- when (f) = and u = 0, the junction is in a 
0-state; by applying a gate voltage u = O.Stt the current 
is strongly suppressed (dotted curve), although its sign is 
always positive; however, when a magnetic flux (j) — O.Stt 
is also introduced, the current sign is reversed (dashed 
curve). Figure [IJb) and Fig.[2Ic) refer to the Josephson 
current for x = 0, and show that its sign can be reversed 
by tuning either the magnetic flux (b) or the gate volt- 
age (c). Notably, the curves of Fig. ^h) demonstrate 
the crucial role of the gate voltage: for m = 0, although 
J switches its sign upon varying the flux, the amplitude 
of the current is very small (solid curve); the latter is en- 
hanced by increasing u (dotted and dashed curves) [3]. 
The curves of Fig. [2fc) refer to different values of mag- 
netic flux; in particular, when no magnetic flux is present 
(solid curve) the TRS yields J(x = 0) = independent 
of the gate voltage u, whereas when 7^ a supercurrent 
can flow and its sign varies with u (dotted curve), with 
a strong enhancement when the flux approaches = tt. 
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FIG. 3: Dimensionless Josephson current vs. u and (f> for an 
AB ring with |A| = 2El, kpL = 60, 7 = 7r/3, and x = 0. 



The whole behavior of the Josephson current J{x = 0) 
as a function of u and 4> is displayed in Fig. [3] The su- 
percurrent exhibits lobes of alternate signs whose nodes 
are located around some special values of cf) and u. The 
latter can be roughly estimated through the condition 
cosxo — —1 (see Eg. p^ l. yielding cj) ~ (2mi + l)7r, 
u ~ (2/712 + l)7r and u ~ 2(7rm3 — kpL), with inte- 
ger m," 's hj. The sign reversal is thus easily controlled 
around these values. This demonstrates the full tunabil- 
ity of the supercurrent through electro-magnetostatic AB 
effect. 

In conclusion, we have proposed Aharonov-Bohm in- 
terferometry as a novel method to realize a controllable 
Josephson 7r-junction. While the magnitude of the su- 
percurrent can be tuned by either the electrostatic or 
the magnetic AB effect, its sign can be controlled only if 
time-reversal symmetry is broken. We have also shown 
that a magnetic field alone does not provide an efficient 
tuning of the supercurrent (see Fig. [2jb)). By contrast, 
the combined use of magnetic and electric field enhances 
the supercurrent amplitude, allowing a full control of the 
junction. In addition, our results also imply that su- 
percurrent measurements can be used to determine the 
asymmetry of the S-matrix; indeed the ordinary DC cur- 
rent in an ring contacted to normal leads only depends 
on the transmission coefficient T — |rp — |t*P, and the 
asymmetry cannot be probed. 
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